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ABSTRACT 

Photoionization cross sections and recombination rate coefficients are presented for the 
L-shell ground state fine structure levels 2s^2p'^ '^P2.o,i of Fe XfX. Several sets of calcu- 
lations including relativistic effects are carried out: (i) Breit-Pauli R-matrix (BPRM), 
(ii) Relativistic Distorted Wave (RDW), and (iii) a semi-relativistic calculation. Non- 
relativistic LS coupling calculations are also done for comparison. The BPRM calcula- 
tions employ a configuration interaction target representation for Fe XX consisting of 
12 LS terms (23-fine structure levels), as in the recently reported BPRM calculations 
by Donnelly et al. (MNRAS, 307, 595, 1999). The background cross sections in all three 
sets of present calculations agree with one another, but differ considerably from those 
of Donnelly et al. . Owing to much more extensive resonance structures in the present 
BPRM calculations, the sum of the corresponding recombination rate coefficients for 
the ^P2,o,i levels are up to 50% higher than the LS rates at low temperarures but 
comparable for higher temperatures; in contrast to the results of Donnelly et al. who 
obtained the LS rates to be higher than their BPRM results by about a factor of 2. 
Reasons for these discrepancies are discussed. 
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1 INTRODUCTION 



Photoionization and recombination are fundamental atomic 
processes in astrophysical plasmas. However, rather elabo- 
rate and extensive atomic calculations are required in order 
to compute the related photoionization and recombination 
coefficients in an accurate and self-consistent manner. In 
most existing astrophysical models these atomic data are 
from methods that do not adequately consider the level of 
complexity attributable to the large number of autoionizing 
resonances that manifest themselves in the photoionization 
cross sections, and consequently, in the inverse process of 
(electron - ion) recombination. Most previous atomic calcu- 
lations divide the recombination process into non-resonant 
'radiative recombination' (RR) and resonant 'di-electronic 
recombination' (DR). This division is valid only if either 
RR or DR is dominant in a given energy region. In gen- 
eral, in the low energy region, the electron-ion recombina- 
tion may contain significant contributions from both the 
autoionizing resonances and the background cross section. 
Therefore, an accurate theoretical treatment requires both 
the RR and DR to be treated in a unified manner. Recent de- 
velopments enable such a unified and self-consistent calcula- 
tion of photoionization and recombination cross sections, in- 
cluding highly detailed resonance structures, using the close 
coupling approximation and the R-matrix method (Nahar 
& Pradhan 1992, 1994; NPl and NP2). For highly charged 



ions, relativistic effects may be considered in the Breit-Pauli 
R-matrix (BPRM) approximation (Zhang & Pradhan 1997; 
Zhang, Nahar & Pradhan (ZNP) 1999). Electron- ion recom- 
bination rate coefficients are thereby obtained including the 
RR and DR processes through ab initio calculations. 

In high temperature X-ray sources K- and L-shell 
atomic ions are prominent constituents. In particular the 
highly ionized ions of iron are important in laboratory and 
astrophysical sources (Liedahl etal. 1992). Recently, Don- 
nelly etal. (1999) employed the BPRM approximation to 
calculate total and partial photoionization cross sections 
for the ground state fine structure levels 2s^2p^ ^^2,0,1 of 
Fe XIX. They also calculated recombination rate coefficients 
for the transitions from the ground level 2s^2p"^ *'S'3/2 in 
Fe XX to those levels in Fe XIX using the partial photoion- 
ization cross sections. They employed a sophisticated repre- 
sentation for the 12 lowest SLtt terms corresponding to 23 
Jn fine structure levels in the target ion Fe XX. Donnelly 
etal. obtained fine structure photoionization cross sections 
for the three ground state levels such that the individual 
cross sections differed considerably from each other, in back- 
ground and in details of resonances, and the weighted aver- 
age differed from the non-relativistic LS values by more than 
a factor of two. The results were surprising since they im- 
ply relativistic efi'ects are larger than expected. Otherwise, 
the individual fine structure cross sections should have a 
background of the same magnitude (although different res- 
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onance structures). For example, an earlier BPRM study 
on C-Iike Fe XXI, with the same ground state symmetries 
^0,1,2 showed that the background cross sections for the 
three levels were nearly equal (Zhang 1998) . We have there- 
fore repeated the calculations of Donnelly et al. as precisely 
as possible, using the same eigenfunction expansion as in 
their work. We obtain very different results, confirmed by 
several other calculations in different approximations, as 
well as much smaller relativistic effects. New recombination 
rate coefficients for the Fe XIX ground levels are also pre- 
sented and compared. 



2 THEORY 

The extension of the close coupling method to electron- 
ion recombination is described in earlier works (NPl, NP2, 
ZNP), together with the details of the unified treatment. 
Here we present a brief description of the theory relevant 
to the present calculations. The calculations are carried out 
in the close coupling (CC) approximation employing the R- 
matrix method, in LS coupling and in intermediate coupling 
with the BP Hamiltonian. The target ion is represented by 
an A^-electron system, and the total wavefunction expan- 
sion, '^{E), of the (A'^-l-1) electron-ion system of symmetry 
SLtt (LS coupling) or Jtt (BP, intermediate coupling) may 
be represented in terms of the target eigenfunctions as; 



*(ij) = A^x.e> + X]' 



,<E>, 



(1) 



where Xi is the target wavefunction in a specific state SiLiiti 
or JiTTi and 9i is the wavefunction for the (A'^+l)-th electron 
in a channel labeled as SiLi{Ji)TViki £i{S Ln or Jiv); kf being 
its incident kinetic energy. $ j 's are the correlation functions 
of the (A'^-l-l)-electron system that account for short range 
correlation and the orthogonality between the continuum 
and the bound orbitals. Most of the present calculations are 
in the BP approximation (Pradhan & Zhang 1997, Zhang 
1998). 

Recombination of an incoming electron to the target ion 
may occur through non-resonant, background continuum, 
usually referred to as radiative recombination (RR): 



e + X^ 



hv + X^ 



(2) 



which is the inverse process of direct photoionization, or 
through the two-step recombination process via autoionizing 
resonances, i.e. dielectronic recombination (DR): 



e + X^ 



(X^ 



e + X++ 
hu + X+ 



(3) 



where the incident electron is in a quasi-bound doubly- 
excited state which leads either to (i) autoionization, a 
radiation-less transition to a lower state of the ion and the 
free electron, or to (ii) radiative stabilization predominantly 
via decay of the ion core, usually to the ground state, and 
the bound electron. 

In the unified treatment the photoionization cross sec- 
tions, CTpi, of a large number of low-n bound states - all 
possible states with n < Wmax ~ 10 - are obtained in 
the close coupling (CC) approximation as in the Opacity 
Project (Seaton 1987). Coupled channel calculations for api 
include both the background and the resonance structures 



(due to the doubly excited autoionizing states) in the cross 
sections. The recombination cross section, ctrc, is related to 
crpi, through detailed balance (Milne relation) as 



a^ 3j (e + Ipf 
4 gi e 



(4) 



in Rydberg units; gi and gj axe the statistical weight fac- 
tors of the recombining (target or residual) and recombined 
ions respectively; a is the fine structure constant, e is the 
photoelectron energy, and 7p is the ionization potential. In 
the present work, it is assumed that the recombining ion is 
in the ground state, and recombination can take place into 
the ground or any of the excited recombined (e-|-ion) states. 
The contributions of these bound states to the total ctrc are 
obtained by summing over the contributions from individ- 
ual cross sections, ctrc thus obtained from ctri, including the 
autoionizing resonances, corresponds to the total (DR-I-RR) 
unified recombination cross section. 

Recombination into the high-n states, i.e. rimax < n < 
ex, must also be included (Fig. 1 of NPl). To each ex- 
cited threshold SiLi{Ji)-Ki of the A'^-electron target ion, 
there corresponds an infinite series of (A''-|-l)-electron states, 
SiLi[Ji)-Kivl, to which recombination can occur, where v is 
the effective quantum number. For these states DR domi- 
nates the recombination process and the background recom- 
bination is negligibly small. The contributions from these 
states are added by calculating the collision strengths, Qdr, 
employing the precise theory of radiation damping by Bell 
& Seaton (1985). Several aspects related to the application 
of the theory to the calculation of DR collision strengths are 
described in the references cited. 

The total (e -I- ion) recombination rate coefficients may 
be calculated using the photoionization cross sections CTpi(e) 
and the DR collision strengths Odr. In the present work 
however, we confine ourselves to a calculation of the level- 
specific ctrc(T, i -^ j), where j refer to the ground state 
fine structure levels ^P2,o,i- The corresponding maxellian 
averaged rate coefficients are obtained as 



a-Rc{T-i -^ j) 



31. 



(5) 



9i kTcV2nm'^kT 
where E = hv = e + Ip is the photon energy. 

3 COMPUTATIONS 

3.1 Residual ion representation for Fe XX 

Close coupling calculations for an (e + ion) system require 
an accurate wavefunction expansion for the states of inter- 
est in the target or the residual ion. Donnelly etal. (1999) 
used an elaborate CI representation for the Fe XX states ob- 
tained using the atomic structure code CIVS (Hibbert 1975). 
We obtain a very similar CI target expansion using the SU- 
PERSTRUCTURE program (Eissner etal. 1974), with the 
same target states as Donnelly etal. , i.e. the ones domi- 
nated by configurations 2s^2p'^,2s2p'^,2p^,2s^2p^3s and all 
associated terms and fine structure levels. Both the CIVS 
and the SUPERSTRUCTURE are carried out in interme- 
diate coupling using the Breit-Pauli approximation. Table 1 
presents the fine-structure level energies from the present 
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Table 1. Calculated (cal) LSJ energy levels of Fe XX in Ry- 
dbcrgs and comparison with the observed (obs) energies from 
NIST (1985), and with those (DBK) by Donnelly etal. (1999). 



Level 



-£^cai(LSJ) 



E{DBK) 



1 


2s22p3 450 


3/2 


0.0000 


0.0000 


0.0000 


2 


2s22p3 ^D° 


3/2 


1.2789 


1.2632 


1.2887 


3 


2s22p3 2£,o 


5/2 


1.6446 


1.6050 


1.6612 


4 


2s'22pS 2po 


1/2 


2.3700 


2.3718 


2.3881 


5 


2s22p3 2po 


3/2 


2.9332 


2.9465 


2.9616 


6 


2s2p4 4p 


5/2 


6.8499 


6.8594 


6.8674 


7 


2s2p4 4p 


3/2 


7.4451 


7.4799 


7.4657 


8 


2s2p4 4p 


1/2 


7.6432 


7.6796 


7.6629 


9 


2s2p4 2£, 


3/2 


9.5480 


9.5006 


9.5729 


10 


2s2p4 2£, 


5/2 


9.6762 


9.6445 


9.7304 


11 


2s2p4 25 


1/2 


10.9215 


10.8920 


10.9652 


12 


2s2p4 2p 


3/2 


11.3899 


11.3219 


11.4352 


13 


2s2p4 2p 


1/2 


12.2595 


12.2113 


12.3052 


14 


2p5 2po 


3/2 


17.8936 


17.8109 


17.9667 


15 


2p5 2po 


1/2 


18.8721 


18.7922 


18.9503 


16 


2s22p23s 4p 


1/2 


65.8494 


- 


65.3888 


17 


2s22p23s 4p 


3/2 


66.3955 


- 


65.9709 


18 


2s22p23s 4p 


5/2 


66.7934 


- 


66.3752 


19 


2s22p23s 2p 


1/2 


66.6709 


- 


66.2395 


20 


2s22p23s 2p 


3/2 


67.0759 


- 


66.4774 


21 


2s22p23s 2d 


5/2 


67.9909 


- 


67.5913 


22 


2s22p23s 2d 


3/2 


68.0771 


- 


67.6739 


23 


2s22p23s 25 


1/2 


69.4502 


- 


68.7533 



BPRM calculation, which are almost identical to those from 
the SUPERSTRUCTURE resuhs, together with experimen- 
tal energies (Sugar & Corliss 1985) and those calculated by 
Donnelly etal. using CIV3. The present calculated energies 
in Table 1 are taken after intermediate coupling stage for 
the target ion in the code RECUPD of the BPRM pack- 
age of codes (Scott and Taylor 1982, Hummer etal. 1993, 
Berrington etal. 1995). As is readily seen, the agreement 
between the calculated and observed energies is very good, 
better than 1% for most of the levels (the present calculated 
energies are slightly better in agreement with experiment 
than those of Donnelly et al. ) Table 2 presents a subset of 
the oscillator strengths calculated, for the transition array 



2s2p , and compared with those of Donnelly 



2s^2p^ 

etal. , and from the U.S. National Institute for Standards 
and Technology (NIST) compilation by Fuhr etal. (1988) 
based on the relativistic calculations by Cheng etal. (1979). 
All three sets of oscillator strengths are in good agreement; 
the present values are in somewhat better agreement with 
the NIST tabulation than Donnelly etal. for some transi- 
tions. We may therefore deduce that the wavefunction rep- 
resentations in the two calculations, Donnelly et al. and ours, 
are of the same accuracy. 



Table 2. Dipole oscillator strengths ( > 0.0001) for fine struc- 
ture transitions in the transition array 2^2 2p'^ -^ 2s2p'^ for 
the target ion Fe XX; comparisons are with the NIST data by 
Cheng etal. (in Fuhr etal. 1988), and with those (DBK) by 
Donnelly etal. (1999). 



Transition 



Present NIST 



DBK 



4 QO 

•^3/2 


- ^Pl/2 


0.0208 


0.0221 


0.0208 




^ ^P3/2 


0.0392 


0.0413 


0.0391 




4p , 
-n5/2 


0.0503 


0.0520 


0.0485 




- ^5i/2 


0.0009 


0.0010 


0.0011 




- '^3/2 


0.0040 


0.0045 


0.0045 




- =^^3/2 


0.0021 


0.0026 


0.0028 


2 no 
-^1/2 


- ^Pr/2 


0.0010 


0.0011 


0.0012 




- ^Si/2 


0.0622 


0.0640 


0.0615 




^ ^Pl/2 


0.0057 


0.0057 


0.0045 




- ^P3/2 


0.0272 


0.0284 


0.0279 




- 'D^/2 


0.0141 


0.0146 


0.0130 


2po 

3/2 


- *^3/2 


0.0010 


0.0011 


0.0011 




- 'P,/2 


0.0003 


0.0003 


0.0003 




- ^Si/2 


0.0026 


0.0030 


0.0018 




- ^Pl/2 


0.0679 


0.0700 


0.0677 




2 p 

-^3/2 


0.0169 


0.0167 


0.0162 




- ^D,/, 


0.0018 


0.0020 


0.0016 




- '^5/2 


0.0242 


0.0250 


0.0236 


2j-)o 


- ^^^1/2 


0.0004 


0.0005 


0.0005 




- "^3/2 


0.0003 


0.0004 


0.0003 




^ 'Ps/2 


0.0033 


0.0038 


0.0041 




- ^5i/2 


0.0289 


0.0300 


0.0297 




- ^Pl/2 


0.0148 


0.0151 


0.0132 




- ^^3/2 


0.0187 


0.0181 


0.0159 




- ^^3/2 


0.0742 


0.0780 


0.0745 


2t-)o 


- "^3/2 


0.0001 


0.0001 


0.0002 




- 'P,/2 


0.0011 


0.0012 


0.0013 




- 'P3/2 


0.0868 


0.0890 


0.0853 




- 2^5/2 


0.0599 


0.0630 


0.0589 



han (1998) and Nahar (1999). Donnelly etal. also provide a 
detailed discussion of the R-matrix method and the calcula- 
tions. In the present work we confine the calculations to the 
ground term levels considered by Donnelly et al. for both the 
total and partial photoionization cross sections. 

In addition to the BPRM calculations using the tar- 
get expansion described in the above sub-section, non- 
relativistic LS coupling calculations, using the same target, 
were carried out for comparison. We also obtained two other 
data sets for the photoionization of Fe XIX, with the rela- 
tivistic distorted wave (RDW) code (Zhang 1998), and with 
the semi-relativistic atomic structure code CATS (Abdallah 
etal. 1988) and ionization code GIPPER (Clark & Abdal- 
lah, unpublished) of the Los Alamos National Laboratory 
(LANL). 



3.2 Photoionization 

The BPRM method as applied to self-consistent photoion- 
ization and recombination calculations has been discussed 
in detail in several previous works and references therein 
(e.g. ZNP, Pradhan & Zhang 1998). The use of such data to 
astrophysical applications in photoionization modeling, and 
photoionization and recombination data for all ions of car- 
bon, nitrogen, and oxygen are described by Nahar & Prad- 



3.3 Recombination rate coefficients 

Maxwellian averaged recombination rate coefficients for the 
ground term levels 2s^2p* ^^2,0,1 are obtained using the par- 
tial photoionization cross sections api for the transitions 
hu + 2s^2p* {^Pj) — > e-|-2s^2p^(*S'3/2). The recombination 
rates were also calculated using the LS coupling photoioniza- 
tion cross sections for hv + 2s'^2p* (^P) — > e + 2s'^2p^{^S°), 
for comparison with the relativistic calculations. 
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3.4 Oscillator strengths 

Another indication of the accuracy of the bound-free pho- 
toionization calculations may be obtained by calculating 
the bound-bound radiative transition probabilities using the 
BPRM method for the fine structure El (dipole allowed 
and intercombination) transitions (e.g. Nahar and Pradhan 
1999a,b). As the same bound state wavefunction expansions 
are employed in both the bound-free and the bound-bound 
transitions, one expects the same accuracy. A number of os- 
cillator strengths of Fe XIX were computed and compared 
with other accurate calculations using the multiconfigura- 
tion Dhac-Fock (MCDF) method by Cheng etal. (1979). 



4 RESULTS AND DISCUSSION 

Results for photoionization, recombination, and transition 
probabilities calculated in the present work are consistent 
with one another, as discussed below. 



4.1 Photoionization cross sections 

Fig. 1 presents the total photoionization cross sections for 
the three ground term fine structure levels ^P2,o,i using the 
BPRM method (the solid hues) and the RDW method (the 
dashed lines). The non-relativistic LS coupling photoioniza- 
tion cross sections of the '^P term and the weighted average 
of the RDW results are also presented (Fig. Id). Compar- 
ing our Fig. 1 with Figs. 1 and 2 in Donnelly etal. one sees 
the present BPRM results differ considerably with those by 
Donnelly et al. in the following three aspects: 1) while the 
results of Donnelly et al. exhibit large variations in the back- 
ground cross sections for the three levels, the present BPRM 
results show those to be nearly the same; 2) in the photon 
energy range ~ 126 — 154 Ry, the cross sections for lev- 
els ^P2 and ^Pi in Donnely etal. 's Fig. 1 show extensive 
and strong resonance structures while our results do not; 3) 
the weighted average (not shown directly here) of the back- 
ground cross sections for the three fine structure levels in 
our results would be very close to the non-relativistic LS 
results, as is obvious from Fig. 1, while those in Donnelly 
et al. differ from the non-relativistic LS values by more than 
a factor of two. 

However, the detailed resonance structures in the 
present BPRM calculations are more extensive and differ sig- 
nificantly from the LS calculations. Fig. 2 shows the cross 
sections at higher resolution. It is apparent that the reso- 
nance positions in the LS "^P results (panel d) do not corre- 
spond in detail or magnitude to the individual resonances in 
the cross sections for the fine structure components '^P2,o,i 
(panels a,b,c respectively). 

Fig. 3 presents the partial cross sections for photoion- 
ization from the ground level ^P2 of Fe XIX into a few fine 
structure levels of the residual ion Fe XX. Comparing this 
figure with Fig. 3 in Donnelly etal. , one again sees that 
similar significant differences exist in both the background 
and the resonances, especially the resonances in the photo- 
electron energy region 20 - 45 Ry. 

Now we discuss the above three major differences be- 
tween the present results and those in Donnelly etal. 1) 
From our many other calculations, such as those shown in 



Fig. 1 for Fe XXII and Figs. (4c) and (4d) for Fe XXIV 
in Zhang (1998), the relativistic effects do not appear to 
affect the background photoionization cross sections very 
much, and the backgrounds for the fine structure levels cor- 
responding to the same energy term appear to be almost 
the same. Those BPRM results, as well as the present re- 
sults, were all confirmed by the RDW calculations (shown 
as dotted or dashed lines in the corresponding figures) . The 
backgrounds in the present BPRM results also agree with 
the semi-relativistic results from the LANL codes (basically 
similar to the RDW results and therefore not shown). 2) The 
resonances in the photon energy region ~ 126 — 154 Ry, or 
the photoelectron energy region ~ 20 — 45 Ry, obtained by 
Donnelly et al. would correspond to the process 



hiy 



2s' 2p* 



2tMn'l - 



2s^2p^ + e. 



(6) 



with n' > 3; it is obvious that the first step in this pro- 
cess involves two-electron jumps and is not allowed (there 
would be some contributions through the configuration mix- 
ing among the n = 2 and 3 states but these would not 
produce such strong resonances) . 3) As is clear from the dis- 
cussion in 1) , the background in the weighted average of the 
results for the three levels should not differ much from the 
non-relativistic LS-coupling results. This should be true for 
ions with low and intermediate nuclear charge Z, where the 
relativistic effect should not be significant enough to affect 
the background. Indeed, we usually use this fact to check 
our calculations. 

The detailed resonance structures on the other hand 
may differ considerably since they correspond to different 
target thresholds in the BPRM and the LS coupling calcu- 
lations. In the relativistic calculations the Rydberg series of 
resonances are more extensive owing to the many more fine 
structure level splittings in the target than the LS terms. 

We also note that a preliminary set of BPRM calcula- 
tion with a smaller number of target terms, excluding the 
2s^2p^3s, also yielded essentially the same results for the 
photoionization cross sections. This is expected for highly 
charged ions since the configuration interaction between dif- 
ferent n-complexes decreases with ion charge. 

4.2 Recombination rate coefficients 

Table 3 presents the QfR(r) calculated using the ctrc 
obtained from the partial photoionization cross sections 
o-pi(^Pj -^ "Sg/a), for the three levels J = 2, 0, 1. The sum 
of the recombination rate coefficients for these three fine 
structure levels is also given, together with the LS coupling 
results. It is important to note that the total recombination 
rate coefficient into a LS term is a simple sum of the val- 
ues for the individual fine structure levels corresponding to 
this term. This is because the expression for the recombi- 
nation rate (see our Eq. (5) or Eq. (7) in Donnelly etal. ) 
already includes a multiplication by the statistical weight 
factor, which is qj = (2 J + 1) in our Eq. (5) and luj in Don- 
nelly etal. 's Eq. (7). We note that, as shown by Eq. (3) 
of NP2, the total recombination rate coefficient for a given 
target level is obtained by sum,ming over contributions from 
all possible states or levels in the recombined ion. Donnelly 
et al. , on the other hand, performed a "weighted average" in 
order to obtain the recombination rate coefficient for the LS 
term - leading to a further error in their numerical values 
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Table 3. The present BPRM recombination rates as func- 
tions of temperature for the transitions from the ground level 

ground term 2s'^2p^ '^Pj 

(SUM) is also shown and compared with the non-relativistic 

LS-coupling results (LS). 



■"S" 2 in Fe XX to the three levels of Fe XIX in its 
2o„4 3 p,. ^[jg sum of the values for these levels 



Table 4. Dipole oscillator strengths for fine structure tran- 
sitions in the transition array 2s^2p'* ^P — > 2s2p^ 3,lpo £qj, 
Fe XIX; comparisons are with the NIST data by Cheng et al. 
(in Fuhr et al. 1988 



r(K) 



'P2 



'Po 



'Pi 



SUM 



LS 



5.0[4] 
6.0[4] 
7.0 [4] 
8.0[4] 
9.0[4] 
1.0[5] 
2.0[5] 
3.0[5] 
4.0[5] 
5.0[5] 
6.0[5] 
7.0[5] 
8.0[5] 
9.0[5] 
1.0[6] 
2.0[6] 
3.0[6] 
4.0[6] 
5.0[6] 
6.0[6] 
7.0 [6] 
8.0[6] 
9.0[6] 
1.0[7] 



2.57 
2.31 
2.12 
1.96 
1.83 
1.71 
1.04 
7.44 
5.81 
4.77 
4.05 
3.52 
3.12 
2.80 
2.55 
1.38 
9.63 
7.43 
6.03 
5.04 
4.30 
3.73 
3.28 
2.91 



11] 
11] 
11] 
11] 
11] 
11] 
11] 
12] 
12] 
12] 
12] 
12] 
12] 
12] 
12] 
12] 
13] 
13] 
13] 
13] 
13] 
13] 
13] 
13] 



4.24 
3.99 
3.77 
3.55 
3.36 
3.18 
2.01 
1.45 
1.13 
9.33 
7.93 
6.90 
6.12 
5.50 
5.01 
2.70 
1.89 
1.45 
1.18 
9.84 
8.40 
7.29 
6.40 
5.69 



12] 

12] 
12] 
12] 
12] 
12] 
12] 
12] 
12] 
13] 
13] 
13] 
13] 
13] 
13] 
13] 
13] 
13] 
13] 
14] 
14] 
14] 
14] 
14] 



1.11 
1.03 
9.63 
9.06 
8.55 
8.09 
5.15 
3.75 
2.94 
2.43 
2.07 
1.80 
1.60 
1.43 
1.30 
6.98 
4.86 
3.73 
3.02 
2.52 
2.15 
1.86 
1.63 
1.45 



11] 
11] 
12] 
12] 
12] 
12] 
12] 
12] 
12] 
12] 
12] 
12] 
12] 
12] 
12] 
13] 
13] 
13] 
13] 
13] 
13] 
13] 
13] 
13] 



4.11 
3.74 
3.46 
3.22 
3.02 
2.84 
1.75 
1.26 
9.88 
8.12 
6.91 
6.01 
5.33 
4.79 
4.35 
2.34 
1.64 
1.26 
1.02 
8.54 
7.29 
6.32 
5.55 
4.93 



■11] 
■11] 
■11] 
■11] 
■11] 
■11] 
■11] 
■11] 
12] 
12] 
12] 
12] 
12] 
12] 
12] 
12] 
12] 
12] 
12] 
13] 
13] 
13] 
13] 
13] 



2.78 
2.57 
2.39 
2.24 
2.10 
1.98 
1.25 
9.15 
7.24 
6.01 
5.16 
4.54 
4.06 
3.68 
3.38 
1.94 
1.41 
1.11 
9.12 
7.70 
6.62 
5.78 
5.10 
4.54 



11] 
11] 
11] 
11] 
11] 
11] 
11] 
12] 
12] 
12] 
12] 
12] 
12] 
12] 
12] 
12] 
12] 
12] 
13] 
13] 
13] 
13] 
13] 
13] 



- and found a difTerent kind of discrepancy compared with 
the LS coupling results. 

The present rate coefficients for the '^P term obtained 
by summing over the fine structure values are also signifi- 
cantly higher than the LS coupling results, especially for the 
low temperatures, the BPRM results being larger by up to 
50%. In contrast, the BPRM results in Donnelly etal. are 
up to a factor of 2 smaller than the LS values. The differ- 
ences between the BPRM and the LS coupling values in the 
present results are due to the different details of the reso- 
nance structures in the two cases. Whereas in LS coupling 
the fine structure thresholds belonging to a given LS term 
are degenerate, these are explicitly included in the BPRM 
calculations at the appropriate energies and consequently 
the resonances are more numerous. As displayed earlier in 
Fig. 2, the resonance positions and shapes are also differ- 
ent in the BPRM and LS coupling calculations. Further, 
more levels of the autoionizing terms which do not autoion- 
ize in the LS coupling case contribute in the BPRM case 
(for example, see the discussion on the Is2p3p autoionizing 
states in Sec. TUB in Zhang (1998)). At higher tempera- 
tures the differences between the present LS coupling and 
BPRM rates are smaller, since the overall background cross 
sections are of similar magnitude at higher energies, above 
all n = 2 target thresholds. It might also be pointed out that 
the autoionizing resonances are unlikely to be significantly 
dampened due to radiative decays, to affect the overall cross 
sections and rates, since they involve A n — transitions 
with much smaller radiative rates than autoionization rates. 



Transition 



Present NIST 



'P2 



'Po - 



Spo 


0.129 


0.147 


Ipo 


0.0313 


0.0355 


''p° 


0.296 


0.34 


■Apo 


0.0772 


0.087 


1 no 

-^1 


0.0037 


0.005 


3 no 


0.0898 


0.1035 


■Apo 


0.0617 


0.0705 


Ipo 


0.0029 


0.003 


'■'P° 


0.0986 


0.1122 



Further discussion of recombination theory in the close cou- 
pling formulation, including the R-matrix method, is given 
in several previous works: Davies and Seaton (1969), Bell 
and Seaton (1985), Robicheaux etal. (1995), Gorczyca and 
Badnell (1997), Pradhan and Zhang (1997), BadneU etal. 
(1998), Robicheaux (1998), Seaton (1998), and Zhang etal. 
(1999). 

4.3 Transition Probabilities 

Transition probabilities provide an independent test of the 
accuracy of the bound state wavefunctions. Fine structure 
transition probabilities may be computed for El dipole 
allowed and intercombination transitions using recent de- 
velopoments in the BPRM method (Hummer et al. 1993, 
Berrington etal. 1998, Nahar and Pradhan, 1999a,b). We 
calculate the oscillator strengths for the fine structure com- 
ponents within several multiplets corresponding to the ini- 
tial ^P symmetry and final ^'^ {S, P, D)° symmetries. Table 4 
presents a subset of these results: the gf- values for the tran- 
sitions 2s^2p* ^P2,o,i — > 2s2p'' ^'^ A°2- The present resuhs 
are compared with the MCDF calculations by Cheng etal. 
(1979) that have been incorporated in the evaluated com- 
pilation by the U.S. National Institute of Standards and 
Technology (NIST; Fuhr etal. 1988). In all cases the agree- 
ment between the present BPRM results and the MCDF 
values is within about 10%. Both the length and the ve- 
locity oscillator strengths are obtained, and found to be in 
good agreement for the stronger transitions (gf > 0.01) in 
all multiplets calculated. 
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Photoionization and recombination of Fe XIX 
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Figure 1. Total photoionization cross sections for tlie three levels 
J = 2, 0, 1 (a, b and c) of the ground term 2s'^2p''' 3p in Fe XIX 
by the present BPRM (solid lines) and the RDW (dashed lines) 
calculations. The non-relativistic LS coupling results are also pre- 
sented (d) , together with the average values of the RDW results 
(the dashed line). 
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Figure 2. As in Fig. 1, but witii resonance structures at liigli- 
resolution. 
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Photoionization and recombination of Fe XIX 9 
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Figure 3. Partial photoionization cross sections for tlie transi- 
tions from the ground level 2s^2p* ^P2 iu Fe XIX to one of the 
low-lying levels in the residual ion Fe XX by the present BPRM 
calculation. 
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